Optical bistability and hysteresis of hybrid metal-semiconductor nano-dimer 
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Optical response of an artificial composite nano-dimer comprising a semiconductor quantum dot 
and a metal nanosphere is analyzed theoretically. We show that internal degrees of freedom of the 
system can manifest bistability and optical hysteresis as functions of the incident field intensity. 
We argue that these effects can be observed for the real world systems, such as a CdSe quantum 
dot and an Au nanoparticle hybrids. These properties can be revealed by measuring the optical 
hysteresis of the Rayleigh scattering. We show also that the total dipole moment of the system 
can be switched abruptly between its two stable states by small changes in the excitation intensity. 
The latter promises various applications in the field of all-optical processing at nanoscale, the most 
underlying of them being the volatile optical memory. 

PACS numbers: 78.67.-n 



I. INTRODUCTION 

Arrays of metallic nanoparticles (often referred to as 
plasmonic arrays), are widely recognized as potential 
building blocks for nanoscale optical circuitry 1 9 (sec also 
Ref.[l3for an overview). Recently, a number of groups re- 
ported fascinating properties of artificial molecules com- 
prised of a semiconductor quantum dot (SQD) in the 
proximity of a metallic nanoparticles (MNP)*^— Non- 
linear Fano resonance s] 12 ' 13 and bistability in the ab- 
sorption spectrum ) 13 ' 1 control of the exciton emission 
of the SQD (inhibition or enhancement)^ and variable 
quenching of the SQD photoluminescence by proximate 
gold nanoparticles as well as "meta-molccular " res- 
onances^ the inhibition of optical excitation and en- 
hancement of Rabi flopping, 1 - tunable nanoswitching^ 
and gain without inversion^ have been predicted. The 
role of the multipolc SQD-MNP interaction in explain- 
ing the spectra of hybrid systems has been discussed in 
details in Ref. HtJ All these effects depend on both ge- 
ometrical parameters and material properties of hybrid 
clusters providing an excellent opportunity of more fine- 
grained control of spectral and dynamical properties of 
nanoscale objects. 

We consider the simplest hybrid nano-cluster compris- 
ing a SQD and a spherical MNP - artificial hybrid di- 
atomic nano- molecule. When such a system is excited op- 
tically, the dipole moment of the optical transition in the 
SQD generates an additional electric field at the MNP, 
which is superposed with the external field. Similarly, 
the induced dipole moment of the MNP generates an ad- 
ditional electric field in the SQD. Thus, the presence of 
the MNP leads to a self action (feedback) of the SQD. 
Together with nonlinearity of the SQD itself, this can 
give rise to a variety of new optical properties. In partic- 
ular, if the coupling between two nanoparticles is strong 
enough, the self-action can result in optical bistability of 



the response. Note that a dimcr comprised of strongly 
coupled two-level molecules can not manifest bistabil- 
ityi2i Thus, a SQD-MNP hetero-dimer is a fascinating 
nanoscopic system exhibiting this feature. 

To demonstrate the feasibility of the bistable opti- 
cal response of hybrid composites, we consider a closely 
spaced CdSe (or CdSe/ZnSe) SQD and an Au nano- 
sphere. We show that for a range of geometrical param- 
eters of the system (SQD and MNP radii and center-to- 
center distance), the optical bistability and hysteresis can 
be observed in it. We argue also that because of the axial 
symmetry of such artificial diatomic molecule, its state 
can be switched not only by traditional change of the 
driving field amplitude, but also by changing the incom- 
ing field polarization with respect to the molecule axis, 
which offers an additional mechanism of control. The 
fact that both the SQD and the MNP can sustain high 
electric fields suggests such possible applications of artifi- 
cial molecules as all-optical switches and optical memory 
cells at nanoscale in the visible; the two stable states of 
the systems have different total dipole moments, provid- 
ing a possibility to store information in this degree of 
freedom. 

The paper is organized as follows. In the next Section 
the model and formalism are described. In Sec. IIIII we 
present the standard steady-state analysis of the bistable 
optical response, which gives rise to the optical hysteresis 
addresses in the Section |TV] We discuss possible applica- 
tions of the predicted effects in Sec. |V] while Section IVII 
summarizes the paper. 



II. FORMALISM 

We assume that the SQD-MNP hybrid molecule is em- 
bedded in a dielectric host with the permittivity Eb and is 
driven by a linearly polarized external electric field with 
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Here, S = diag(— 1, — 1,2) is the angular part of the 
dipole field Green's tensor (z axis being parallel to the 
system axis), d is the SQD-MNP center-to-center dis- 
tance while Pmnp is given by: 
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FIG. 1: Schematics of the hybrid SQD-MNP system embed- 
ded into a homogeneous dielectric background with the per- 
mittivity Eb and subjected to an external field with an ampli- 
tude Eo. The SQD optical transition dipole moment and the 
semiconductor dielectric constant are denoted as p, and e s , re- 
spectively, while a(w) is the MNP polarizability. The curved 
arrows symbolize the dipole-dipole SQD-MNP interaction. 



an amplitude E and frequency u. Figure Q] shows the 
schematics of the system. The SQD is modeled as a two 
level system with the transition frequency ujq and opti- 
cal transition dipole moment fi. It is treated quantum 
mechanically within the framework of the Maxwell-Bloch 
equations for the 2x2 density matrix p mn (m, n — 0, 1). 
The MNP is considered classically; the response of the 
MNP is described by its frequency dependent scalar po- 
larizability within the point dipole approximation (this 
can easily be generalized for the case of more complex 
shapes of the MNP by considering an appropriate polar- 
izability tensor). All sizes of the system (the SQD and 
MNP radii and the SQD-MNP ccnter-to-center distance) 
are assumed to be small enough to neglect the retardation 
effects and to consider both particles as point dipolcs. 
The rotating wave approximation is used throughout the 
paper, so that all time dependent quantities represent 
amplitudes of the corresponding characteristics of SQD, 
the set of equations for which reads: 
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where Z = pn — poo is the population difference between 
the SQD excited and ground states, R is the amplitude of 
the off-diagonal density matrix element defined through 
pio = — (i/2)i?cxp(— iuxt), 7 and T are the relaxation 
constant of the population and the dipole dephasing, re- 
spectively, A = luq — uj is the detuning of the driving field 
from the SQD resonance, and SI = fiE/H is the electric 
field (in frequency units) acting inside the SQD. The field 
acting upon the SQD is equal to the sum of the exter- 
nal field E and the field produced by the induced dipole 
moment Pmnp of the MNP. The field E inside the SQD 
is reduced by the factor e' s = (e s + 2e b )/(3e b ) where e « is 
the permeability of the SQD (see e. g. Ref . l22l or Ref. 1231 
Ch. V, p. 138): 
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where a(ui) = a 3 is the classical frequency de- 

pendent polarizability of the MNP, a being its radius, 
7(w) = [em(w) - e b ]/[e M {uj) + 2e b ], and e M {uj) is the 
dielectric function of the metal. We do not take into 
account the corrections to the polarizability due to the 
depolarization shift and radiative damping j2a which are 
negligible for nanoparticle sizes of our interest (<10 ran). 
The second term in the parenthesis of Eq. ([3]) is the field 
produced by the SQD dipole moment Psqd = —ipbR at 
the MNP. 

The exciton radius in CdSc is about 5 nrn^ while the 
typical radius of the considered SQD is about 1.5 — 2 nm, 
so the wavefunctions involved in the optical transition 
are extended over the whole dot. In deriving Eq. (|3J) 
we used therefore the approximation of the homogeneous 
electric polarization of the whole SQD volume. In this 
case the dipole field around the SQD is screened by the 
bare background dielectric constant only^ Note that the 
dipole moment Psqd is calculated quantum mechanically 
and accounts for the screening which results from the 
SQD dielectric response (see below). Finally, for the total 
electric field inside the SQD we obtain: 
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As is seen from Eq. © , the presence of the MNP results 
in two effects: the first term accounts for the renormaliza- 
tion of the external field amplitude Eo while the second 
represents the self-action of the SQD via the MNP; the 
field inside the SQD depends on the dipole moment of 
the SQD itself. 

The effect of the self-action on the dynamics of the 
SQD-MNP hybrid nano-molecule can be revealed after 
substituting Eq. ((4]) into Eq. (jlb[) and representing SI in 
the form 



SI = Sl - i GR 
with Qq and G given by 
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where Qo = A*Eo /H is the Rabi frequency of the bare ex- 
ternal field, Qq is the renormalizcd Rabi frequency and G 



3 



is the feedback parameter. The latter absorbs all infor- 
mation governing the SQD self-action, such as, material 
constants, geometry of the system or details of the inter- 
action (e. g. contributions of higher multipoles^). 

In a number of recent publications dealing with the 
same system, a different formula for the constant G was 
used in which the factor e' s appears squared in the denom- 
inator of the GJ^t— The second factor e' s is supposed to 
take into account the screening of the dipole field by the 
SQD dielectric response. We note, however, that it is the 
product G R that determines the latter field and, as can 
be seen from Eq. ([8b]) . R oc On oc Oo/£g, so the dipole 
field G R is already additionally screened. 

The feedback G is the most important parameter of the 
theory; once it is calculated it determines the nonlinear 
properties of the SQD response. Using Eq. J5]), Eq. JTbJ 
can be rewritten in the following form: 

R = -[(T-GiZ) + i(A + G R Z)}R + n Z (7) 

with G R = Re(G) and G\ = Im(G). From Eq. © two 
consequences of the SQD self-action become apparent: 
(i) - the renormalization of the SQD resonance frequency 
ujq l— ^ + Gr Z and (ii) - the renormalization of the 
dipole dephasing rate r M> Y — G\ Z . Both renormalized 
quantities depend on the population difference Z . Simi- 
lar renormalizations originate from the local field correc- 
tion in the nonlinear optical response of dense gaseous 
assemblies of two-level systems^ optically dense thin 
films,— and linear molecular aggregates.— The popula- 
tion dependencies of the SQD resonance frequency and 
the dipole dephasing rate provide a feedback mechanism, 
which results in a number of fascinating effects. 
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FIG. 2: Two possible types of the SQD excitation. Left plot 
- the excitation frequency uj lies below the renormalized SQD 
transition frequency cjn — Gr (A > Gr); as Z increases with 
the excitation the system is driven further out of resonance. 
Right plot: uj > loo — Gr (A < Gr); the excitation drives 
the system into the resonance which favors the bistability to 
occur. 

Let us assume that both Gr > and G\ > 0, which is 
the case for a CdSc SQD conjugated with a golden MNP. 
Then the renormalized resonance frequency increases 
with the excitation intensity, ranging from ujq — Gr 
(Z = — 1 in the ground state) to uj + Gr (Z = 1 in the 
excited state). In this case the response of the system de- 
pends on the relative position of the excitation frequency 
uj with respect to the renormalized SQD transition fre- 
quency ujq — Gr. Thus, if uj < ujq — Gr or equivalently 



A > Gr (see Fig. [2 left plot) then the excitation is driv- 
ing the SQD out of resonance, so that the SQD is becom- 
ing less absorptive. Contrary to that, if uj > ujq — Gr or 
A < Gr (see Fig. [3J right plot) the SQD is being driven 
into the self-sustaining resonance by the incoming field. 
That is the case of the positive loopback. In this case, 
apart from the usual linear "weak field" solution, the 
second kind of a stable state can turn up, which results 
from the above-mentioned positive feedback mechanism. 
We show below that the latter has a threshold charac- 
ter, giving rise to bistability and hysteresis of the system 
response characteristics. 



III. STEADY-STATE ANALYSIS 



First, we analyze Eq. ([Taj) and Eq. (jlbft under steady- 
state conditions (Z = R = 0) to obtain stationary states 
of the system. The corresponding solutions read: 
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Equation (|8aj) is of the third order in Z and therefore 
may have three real solutions, depending on values of A, 
r, Gr , and G\. The same applies to the SQD dipole 
moment amplitude R. 

Hereafter, wc consider a CdSc SQD in the vicinity 
of an Au MNP and use the following set of parame- 
ters: the transition energy Hujo = 2.36 cV (which cor- 
responds to the optical transition in a 3.3 nm SQD), 
the SQD dielectric constant e s = 6.2, the SQD tran- 
sition dipole moment /j, = 0.65 e-nm,— the MNP ra- 
dius a = 10 nm, the SQD-MNP ccnter-to-ccntcr distance 
d = 17 nm, the host dielectric constant £}, = 1, and the 
relaxation constants 7 and T defined through 1/7 = 0.8 
ns and 1/r = 0.3 nsJ^ To calculate the polarizability 
7(0;) of the MNP, we used the tabulated data for the 
permittivity of gold from Ref. H^. For these parameters 
G = G R +1G1 = (25.4+ 10.6 i)T. Note that the fre- 
quency domain of our interest is a narrow region in the 
vicinity of the SQD resonance, with the width of about 
several units of T (sec below) , which is much smaller that 
the width of the MNP plasmonic resonance. We therefore 
neglected the frequency dependence of the MNP polar- 
izability when calculating the feedback parameter G and 
used 7(w) w 7(^0). 

Figure [3] shows the solution of Eq. ([5a]) and Eq. (|8b)l 
for the set of parameters specified above and different 
detunings A. As is seen from the plots, within a win- 
dow of -3.3 T < A < 13.8 T, the field dependence of Z 
and R, have three allowed values for a given intensity 
No| 2 /(7 r )- Tnc u PP cr limit of the window, A = 13.8T, 
corresponds to the excitation frequency uj = ujq — 13.8 T 
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FIG. 3: Stationary solutions to Eq. (|8al) and Eq. (|8b|) for a 
hybrid nano-molecule comprising a 3.3 nm CdSe SQD and 
20 nm Au sphere separated by the center-to-center distance 
of 17 nm (other parameters are described in the text). The 
solutions are calculated for different detuning A as functions 
of the normalized external field intensity Qo / VtF- The upper 
row shows the population difference Z = pu — poo (left) and 
the SQD dipole moment amplitude \R\ (right). The lower left 
panel displays the absorptive part of the dipole amplitude, 
Im(i?), while the right one — its dispersive part, Re(i?). 



which lies above the renormalized SQD resonance fre- 
quency uq — Gr, = ujq — 25.4 r (the positive loopback 
case). The lower limit of the window, A = — 3.3T, is 
negative and the corresponding frequency lies above the 
bare resonance. Nevertheless, the SQD can still be driven 
into the self-sustaining resonance by the external field. 
At larger u (larger negative A), the resonance between 
the excitation and the SQD can not be attained because 
it requires a significant positive population difference Z 
that is unreachable under stationary conditions: due to 
the saturation effect, the upper limit for the population 
difference is Z = in the steady state. Because of that, 
the bistability effect disappears for large negative detun- 
ings. 

Optical bistability can therefore be observed within a 
window of detunings in the vicinity of the SQD transition 
frequency. The width of the window is typically on the 
order of several units/few dozens of T. In small SQDs 
the dipole dephasing time is usually strongly temperature 
dependent and can change from ns for low temperatures 
to less than ps for room temperature^ The bistability 
exists only if the feedback parameter G exceeds some 
threshold value. If Gj = the condition for bistability to 
occur is G R > 4P2S, while if Gr = the condition is Gi > 
8r£L The value of G is determined by the geometry and 
material properties and can not be increased arbitrarily. 
Therefore, these criteria can be used to choose suitable 



materials and system configuration. 

We note that only the population difference Z mani- 
fests the standard S-shape curves while all the quantities 
related to the SQD dipole amplitude R exhibit more ex- 
otic coiled curves. As we show below, the latter leads 
to completely different types of hysteresis loops for these 
quantities. 
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FIG. 4: Optical hysteresis loops of the population difference 
Z = pu — poo and the SQD dipole moment amplitude R ob- 
tained by solving Eqs. (|la[) and (|lb|) with the external field 
intensity | S7o | 2 /-^/-yT being adiabatically swept back and forth 
across the bistability region (the sweeping speed of the renor- 
malized Rabi frequency Qo / VoT is 2.1 x lCT 3 r) The detuning 
from the resonance is A = 1.5r. All other parameters are as 
in Fig. E 



IV. OPTICAL HYSTERESIS 

We performed time-domain calculations with the ex- 
ternal field intensity \£lo\ 2 / (fyF) being adiabatically swept 
back and forth across the bistability region (the sweep- 
ing speeds are given in figure captions), monitoring the 
evolution of the system to figure out which steady state 
branches are stable. The results are shown in Fig. 2) 
For the population difference Z and the absorptive part 
of the SQD dipole moment Im(i?) we find the standard 
behavior. Upon increasing the applied intensity, the sys- 
tem follows the lower (stable) branch until the intensity 
reaches the critical value at which the system switches to 
the upper branch (which is also stable). Upon sweeping 
the intensity back, the system stays on the upper branch 
and then switches back down to the lower one at the other 
critical intensity, completing the hysteresis loop. The in- 
termediate branch can not be revealed by the adiabatic 
sweeping of the field because it is unstable which can also 
be checked by the standard stability analysis. 
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Both the and the Im(i?) behave in a very differ- 
ent manner, manifesting hysteresis loops with kinks. In 
the case of the \R\ the upper branch is unstable and the 
hysteresis loop is triangular, while the Im(i?) with its un- 
stable lower branch has even more complicated bow-tie 
hysteresis curve. To the best of our knowledge no such 
optical hysteresis loops have neither been predicted nor 
observed so far. 

These unusual properties are also expected to mani- 
fest themselves in the MNP dipole moment because it 
depends on that of the SQD [see Eq. ([3])] and can there- 
fore be switched abruptly as well. 
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(see Fig. |4]), which suggests that the major contribution 
to the scattering is coming from the SQD. To confirm 
the latter we calculated the relative contribution of the 
SQD dipole moment into the total scattering amplitude, 
-Psqd| 2 /|-Psqd + -Pmnp| 2 which turned up to be on the 
order of unity across the whole hysteresis region. Within 
this region the contribution of the MNP to the scattered 
field is typically an order of magnitude less than that of 
the SQD for the chosen set of parameters. 

It is to be noticed that the ratio of the scattering inten- 
sity in the two stable states can be as high as about 25 for 
A = 0. Such high contrast suggests that the system can 
be used as an optical memory cell: lower intensity can 
represent a logical zero while the higher - a logical one. 
The state of the cell can be switched by sweeping the 
excitation power across the bistability region. Another 
possibility to switch the cell is to maintain the power con- 
stant while changing the incident field polarization. The 
latter mechanism is characteristic for this particular sys- 
tem due to its axial symmetry. We point out that such 
memory is a volatile one as it requires constant pumping. 
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FIG. 5: The Rayleigh scattering intensity |Psqd + -Pmnp| 2 in 
units of \fi\ 2 for different values of A (indicated in the figures) 
as a function of excitation intensity. The corresponding power 
density is shown in the upper axis. Sweeping speeds of the 
renormalized Rabi frequency f2o / x/tT in units of F is on the 
order of 10 -3 for all values of A. All other parameters are as 
in Fig. El 



V. RAYLEIGH SCATTERING AND OPTICAL 
STORAGE 

In experiment, the intensity of the Rayleigh scattering 
can be measured. The amplitude of this characteristic is 
known to be proportional to the squared absolute value 
of the total system dipole moment |Psqd + -Pmnp| 2 and 
is therefore expected to manifest bistability as well. We 
plot the amplitude of the Rayleigh scattering in Fig. [5] 
for different values of the detuning A. The figure shows 
that various types of hysteresis curves can be observed: 
the standard loop (for A = 0) as well as more exotic 
triangular loops (A = 1.5 T and A = 8T) and the bow- 
tie one (A = 4.5 T). The two latter types of hystere- 
sis curves are characteristic for the SQD dipole moment 



VI. SUMMARY 

We investigated theoretically the optical response of a 
hybrid "artificial molecule" comprised of a closely spaced 
spherical semiconductor quantum dot (modeled as a two- 
level system) and a metal nanosphere (considered classi- 
cally) , which are coupled by the dipole-dipole interaction. 
The interaction results in self-action of the SQD via the 
MNP, leading to the population dependence of the SQD 
transition frequency and relaxation constant of the SQD 
dipole moment. This provides a feedback mechanism re- 
sulting in several fascinating effects. Thus, we found that 
the system can manifest bistability and optical hystere- 
sis. In particular, the total dipole moment of the system 
can be switched between its two stable states by the in- 
coming field. The latter suggests such possible applica- 
tions as optical memory cells and all-optical switches at 
nanoscalc in the visible. 

Because the SQD-MNP dipole-dipole interaction de- 
pends orientation of the dipole moments of the two par- 
ticles, the switching can be achieved not only by the tra- 
ditional control by the incident amplitude, but also by 
changing the polarization of the incoming field with re- 
spect to the system axis. Our calculations performed for 
typical system parameters, such as those of a CdSe or 
CdSc/ZnSc quantum dot and an Au nanoparticle com- 
plexes, predict the optical bistability of a SQD-MNP ar- 
tificial molecule. The Rayleigh scattering and modern 
methods of single molecule/particle spectroscop y 30 ' 32 " — 
could probably be used to discover the predicted effects 
experimentally. 

To conclude, we considered the simplest diatomic hy- 
brid artificial nano-molcculc. We expect, however, that 
more complicated clusters (such as an SQD surrounded 
by several MNPs, as considered in Ref. Il5l ) can also ex- 
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hibit these effects because in such systems nanoparticlcs 
are just playing the role of a "resonator" and provide 
feedback to the nonlinear two-level system. Anisotropy 
of nanoparticles can also easily be accounted for by using 
an appropriate tensor instead of the scalar polarizability. 
Finally, we note that a very interesting aspect of this kind 
of systems is the direction of the total induced dipole mo- 
ment, which can also be bistable. 
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